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Abstract. We derive an analytical form of the Schmidt modes of spontaneous parametric down-conversion 
(SPDC) biphotons in both cartesian and polar coordinates. We show that these correspond to Hermite- 
Gauss (HG) or Laguerre-Gauss (LG) modes only for a specific value of their width, and we show how 
such value depends on the experimental parameters. The Schmidt modes that we explicitly derive allow 
one to set up an optimised projection basis that maximises the mutual information gained from a joint 
measurement. The possibility of doing so with LG modes makes it possible to take advantage of the 
properties of orbital angular momentum eigenmodes. We derive a general entropic entanglement measure 
using the Renyi entropy as a function of the Schmidt number, K, and then retrieve the von Neumann 
entropy, S. Using the relation between S and K we show that, for highly entangled states, a non-ideal 
measurement basis does not degrade the number of shared bits by a large extent. More specifically, given 
a non-ideal measurement which corresponds to the loss of a fraction of the total number of modes, we can 
quantify the experimental parameters needed to generate an entangled SPDC state with a sufficiently high 
dimensionality to retain any given fraction of shared bits. 



1 Introduction 

Entangled states are one of the most striking predictions of 
quantum mechanics. Their generation and measurement 
not only allows for tests of the theory but may also re- 
sult in novel technological advances. Entangled states form 
the basis of quantum information and quantum comput- 
ing and are essential for applications such as quantum 
cryptography where they are used for distributing secret 
cryptographic keys between two parties |ll2l3j . 

Quantum systems can be entangled in various degrees 
of freedom, e.g. polarisation, conjugate variables such as 
time and energy, position and momentum, and also in 
multiple variables simultaneously (known as hyperentan- 
glement) |4|5| . Recently there has been a lot of interest in 
demonstrating entanglement between spatial modes car- 
rying orbital angular momentum (0AM) [6 . These are 
of particular interest for quantum information protocols, 
such as quantum key distribution (QKD), as they reside 
within an infinite-dimensional, discrete Hilbert space, thus 
allowing large amounts of information to be impressed 
onto a single photon and multiple bits of secret key to 
be extracted for each entangled state. 

Experimentally, pairs of photons entangled in their 
0AM can be reliably produced using spontaneous para- 
metric down-conversion (SPDC) [7] and there has been 
a lot of effort to calculate the precise form of the down- 
converted photons, and hence quantify the ensuing degree 



of entanglement, by calculating the spiral bandwidth of the 
system [8 .9 lOi. Alternatively, the amount of entanglement 
can be quantified by performing a Schmidt decomposition 
of the down-converted state [11^. This gives information 
both about the pairing of the photons and about the de- 
gree of entanglement of the state via the Schmidt number, 
K, which is defined as the average number of non-zero 
coefficients in the Schmidt decomposition. For an entan- 
gled bipartite pure state the maximum possible correla- 
tion will occur when the measurements are performed in 
states which correspond to a Schmidt decomposition of 
the bipartite state [12114115) . 

Modes carrying 0AM are most commonly described in 
polar coordinates by superpositions of Laguerre- Gaussian 
(LG) modes but may, equivalently, be described in carte- 
sian coordinates by superpositions of Hermite- Gaussian 
(HG) modes. In this paper we calculate explicitly the 
Schmidt bases for both coordinate systems of the SPDC 
state under the assumption of gaussian phase matching. 
Such bases can be experimentally implemented for joint 
measurements. We use our results to calculate the Schmidt 
number of the state and its Renyi entropy. This not only 
gives an alternative measure of the entanglement, but also 
allows us, as a special case, to calculate the von Neumann 
entropy of these reduced states, and hence the maximum 
information that can be shared as a function of the exper- 
imental parameters. 
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The paper is organised as follows. In section 2 we de- 
scribe the biphoton state produced using SPDC. We show 
its dependence on the pump parameters and crystal length 
and describe the gaussian phase matching approximation 
that can be made in the near coUinear limit. This allows 
us, in section 3, to calculate the Schmidt decomposition 
in both cartesian and polar coordinates, and to show how 
one can easily transform between the two. In section 4 we 
use this result to analyse the entanglement by means of 
the Schmidt number and hence the Renyi and von Neu- 
mann entropies. We then demonstrate the simple relation 
between the information entropy and the Schmidt number 
and show how the information loss due to non-ideal mea- 
surements depends on the number of entangled modes and 
how high dimensional entanglement protects the shared 
information. 



If we scale the wavevectors qi.^ by a factor F' = yjhja 
s.t. q = c( / r' then we can re-write ^ as 



V'(q'»,q'J oc exp 



where 



|qi + q ' 

ha 



I 12 



exp (-6cT|q', - q'j2) 



ha ^ 



(4) 



(5) 



with Zj, the Rayleigh range of the pump beam. Writing 
([2]) in this form allows us to see the symmetry between 
the two parts of the wave function and will help explain 
some of our later results. In particular, we can see that 
our results will depend on the product ha and not h and 
a independently. 



2 The SPDC state 

In spontaneous parametric down-conversion (SPDC) two 
lower-frequency photons, commonly referred to as signal 
and idler, are generated when a pump field interacts with 
a nonlinear crystal (TB] . The spatial structure of the down- 
converted biphotons depends both on the pump field and 
on the phase matching. For a gaussian pump it can be 
written [T7]: 

V'(q^, q.) = AAexp |^-^|q. + q»l'^ sine (^^^^ + ^ 

(1) 

where qs,i are the transverse components of the wave vec- 
tors ks_i for the signal and idler fields, Akz is the longitu- 
dinal component of the wave mismatch between the pump 
wave vector kp and the down-converted photons wave vec- 
tors: Akz = (kp - ki - ks)^ ~ (|qi - q^p) /kp, where we 
have made use of the paraxial approximation. <P is an ad- 
ditional phase mismatch which depends on the internal 
refractive indices (i.e. its value can be tuned by tilting the 
crystal or by changing its temperature). For our analysis 
^> = 0. 

Near the coUinear phase matching regime the analysis 
can be simplified by using a gaussian approximation of the 
phase matching term |18ll7j so the state can be written 
in the form: 

V'(q., q.) (X exp j exp {-h^\ci, - q,p) (2) 

where h and a depend on the pump waist, Wp, and wave 
number, kp, and on the crystal length, L, in the following 
way 




3 Schmidt decomposition 

As we stated earlier, a useful, and experimentally conve- 
nient, measure of the entanglement is the Schmidt num- 
ber, K. In order to calculate this we first need to calculate 
the Schmidt decomposition of the down-converted state, 
which is done by writing it in the form 

ai)A\Pi)B, (6) 

i 

where |ai)A, \Pi)B are the Schmidt modes, defined by the 
eigenvectors of the reduced density matrices, and the real 
and positive Schmidt coefficients, are the correspond- 
ing eigenvalues, with each of the factors in the normalized 
set {Aj} representing the probability of detecting the en- 
tangled state in the i^^ entangled Schmidt mode \ ai) A\fii) b ■ 
If all of the coefficients are different, then the Schmidt 
decomposition is unique. Whenever some of the coeffi- 
cients are equal, one has a choice of infinitely many 
different Schmidt bases. 

The Schmidt decomposition provides insights into the 
nature of the bipartite entanglement by determining the 
natural set of biorthogonal mode pairs (or orthonormal 
bases) for the two systems [11I18I19I20I21] while the coef- 
ficients allow us to calculate the Schmidt number, K (i.e. 
the average number of modes in the state) and the entropy 
of entanglement. By knowing the Schmidt decomposition 
explicitly in HG and LG modes, one can easily implement 
an ideal detection basis. 

In this paper we are interested in the biphoton state 
produced using SPDC which is entangled in its 0AM. As 
this can be equivalently described using either HG or LG 
modes we perform the Schmidt decomposition in the two 
corresponding coordinate systems. By approximating the 
biphoton state as the double Gaussian in ^ we are able to 
calculate the Schmidt decomposition in an analytical form 
for both. We also demonstrate their equivalence by using 
the well-known relationship between HG and LG modes 
[22] . We then show that, as expected, the entropic analysis 
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leads to the same conclusions for both. Note that in both 
cases, the Schmidt modes for the SPDC state will have 
the same form for both signal and idler, due to symmetry 
requirements [20j and so we will obtain a decomposition 
of the form: 



(7) 



where the functions Ua,b depend on the coordinate system 
employed, and the labels a and b correspond to different 
degrees of freedom: in the cartesian case a and b will be 
replaced with m and n, in the polar case with £ and p. 
The m and n quantum numbers label the two transverse 
degrees of freedom, while the i and p quantum numbers 
label the angular and radial degrees of freedom, respec- 
tively. The sum is calculated on two indices because we 
perform the decomposition in the two-dimensional plane 
perpendicular to the direction of emission. 



3.1 Decomposition in cartesian coordinates 

In a recent paper, Straupe et al. [23^ reported a proof-of- 
principle experiment demonstrating that an appropriately 
chosen set of HG modes constitutes a Schmidt decomposi- 
tion for transverse momentum states of biphotons gener- 
ated by SPDC. For clarity and completeness we perform 
an equivalent Schmidt decomposition in cartesian coordi- 
nates (details of the calculations are given in Appendix [A| 
before extending our analysis to polar coordinates in tlie 
next section and then demonstrating their equivalence. 

The cartesian decomposition requires a separation of 
each of the variables and qs into a pair of orthogonal 
variables, q and q±, so that the wave function assumes the 
form '0(qi,qs) tp{qi,qs,qi±,qs±)- We define the carte- 
sian basis of HG modes as 



h„{rq) 



(n!2» 0^)1/2 



(8) 



where F = 2 



: Wp Y 2^ is the width of the HG modes. 

If we express the wave function ([2| in terms of this basis 
we obtain 



m,n>0 



where 
and 



hmn{x,v) = h,n{rx)hn{ry) 
b(T- 1 



(10) 



ba + l 



(11) 



Expression ^ is the cartesian Schmidt form for the 
SPDC state with corresponding coefficients: 



(1 - 



'2 \ ,,m+n 



Aba 



{1 + baY 



ba-l 



b(T + 1 



'm-{-n 



(12) 



Note that these are exactly equivalent to those given in 
equation (5) of [23] . 



3.2 Decomposition in polar coordinates 

As LG modes are currently the preferred basis for many 
spatial entanglement experiments we also calculate the 
Schmidt decomposition in polar coordinates (details of the 
calculations are given in Appendix |b|. The polar decom- 
position requires a separation of each of the variables q^ 
and q^ into a pair of polar variables, so that the wave 
function assumes the form ■0(qi,qs) — 'ip{pi, ps,ipi,(ps), 
where p and f are the radial and angular variables. We 
take the LG modes to have the standard definition in mo- 
mentum space, and we add the same scaling factor 7^ as 
in the cartesian case: 



(13) 



where Lp are generalized Laguerre polynomials. One can 
express the wave function ^ in terms of LG modes, of 
width r/V2 



GO OO 



v = (i-m2) J2 E/^''^'"^^p(^'^-^')^Gp'(^^-^^)- 

i— — oo p—0 

(14) 

where T = 2./^ = Wr>\f^. 

The above expression is the polar coordinate form for 
the Schmidt decomposition of the SPDC state, where jjL 
is defined in (11) and the polar Schmidt coefficients are 
given by 



/A^=(l-^2)^2p+K|^ 



Aba 



ba-1 
ba + l 



2p+\i\ 



(15) 



to be compared to ( 12 ). 



3.3 Equivalence of Schmidt bases 

The expressions calculated above are equivalent descrip- 
tions of the entangled state and, just as it is possible to 
transform LG modes into HG modes |24j and vice versa 
[22] , we are also able to convert between our two Schmidt 
bases. In fact, it is straightforward to convert the Schmidt 
decomposition in cartesian coordinates, equation (|8]), into 
that in polar coordinates, (13). The first step is to notice 
that the values of m and n that satisfy m + n — N yield 
the same Schmidt coefficient (1 — p?)p^ \ for the polar 
case this happens for all the values of £ and p that satisfy 
\l\ + 2p ^ N. The number N is called the mode order. 
The relation between the cartesian basis (|8| and the 



polar basis (13) is 



LGl^\p) 



N 

^b'p^k^N-k.kiq^q^) 

fc=0 



(16) 
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where the relation between cartesian and polar coordi- 
nates is the canonical one, and 



AN) 



(17) 



where m + n— \£\+2p = N . This fact enforces the conser- 
vation of the mode order when changing Schmidt basis. 



4.2 Renyi entropy 

An alternative approach to quantifying the entanglement 
of an SPDC state is to calculate its entropy. The most 
famous entropic function is the Shannon entropy, which 
appears in information theory and statistical mechanics 
[25|26l . A more general measure, however, is the Renyi en- 
tropy, which is obtained by neglecting the grouping prop- 
erty of entropy [27128) . For a probability distribution {pk} 
the Reyni entropy of order a is defined as 



4 Analysis of the entanglement 



4.1 Schmidt number 



The entanglement of a state can be quantified by the prob- 
ability distribution of the modes it contains. Intuitively, a 
state is more entangled whenever this probability distribu- 
tion is more 'spread out'. A particularly important mea- 
sure of entanglement is the Schmidt number, K , which 
corresponds to the number of significant modes in the 
Schmidt decomposition |20I21| . For states in the form Q , 
this is defined as 



K 



Tr [p] ■ 
Tr[p2 



1 



(18) 



where p is the reduced state formed by tracing over one 
part of a pure bipartite state and yjXa,b are the Schmidt 
coefficients that appear in the Schmidt decomposition of 
the bipartite state. One can immediately see that the 
eigenvalues of the reduced state are just the square of the 
Schmidt coefficients. A state will be separable when K ~1 
and entangled if i^T > 1. Applying this to the Schmidt de- 
compositions calculated earlier, gives 



K = 



oo oo 
_m— n— 



n -1 



2 



oo oo 



i— — oo p— 



4 V ba 



(19) 



This result agrees with previous calculations of K [T8^ and 
is independent of the Schmidt basis used. This means that 
one has the freedom to choose the basis that best matches 
the experimental conditions with no consequence on the 
dimensionality of the Hilbert space that is spanned by the 
detection basis. 

The effect of the experimental parameters can be seen 
more clearly if, as in ([H]), we write ha = Wp^/Ljkp = 

yjL/2zr where L is the crystal thickness, Wp and kp are 
the width and wave number, respectively, of the pump 
and Zr is its Rayleigh range. Note that K — \, which 
means that the state is not entangled, whenever 6cr = 1, 
which corresponds to choosing experimental parameters 
such that the crystal length is twice the Rayleigh range 

{L = 2Zr). 



Ho.{{Pk}) = 




a > 0. 



(20) 



Note that when a — )■ 1 one regains the Shannon entropy. 

A simple calculation shows that the Renyi entropy of 
the Schmidt coefficients in either eq. ( 12 1 or ( 15 ) is 



— logn ■; 

1 ^2 1 



ba\ 



\\-ha\ 



(21) 



Using equation ( 19 1 we can rewrite the parameter ha in 



terms of the Schmidt number, K,asba = \/K — \/K— 1 
and hence find the Renyi entropy in terms of the Schmidt 
number. Replacing ba with \/ K — \JK —\ in (21) gives 
the Renyi entropy as a function of if, which can be ap- 
proximated by 



iJ„(X)~log2X-/(a) 



(22) 



where /(a) = 2 — log2(a2)/(a — 1). For a discussion of 
the quality of this approximation, see the next subsec- 
tion. This tells us that, to a good approximation, valid for 
sufficiently high K , different orders of the Renyi entropy 
differ by a constant value. Note that /(2) = and thus 
Hi ~ log2 X, in fact the Schmidt number is related to the 
Renyi entropy of order 2 by = 2^^ . 



4.3 Von Neumann entropy 

For quantum systems we can write the quantum Renyi 
entropy as 



1 



1 - a 



log2[Tr(p")], a>0. 



(23) 



where p is the reduced density matrix |29j . An important 
special case of the Renyi entropy is when one takes the 
limit a — 1 in which case ([23]) reduces to the von Neu- 
mann entropy of the reduced state |llll9j 



S'[/5] = -Tr[/51og2p]. 



(24) 



The entropy of a reduced state is known as either the index 
of correlation [12j or the entanglement entropy |30)31j and 
the importance of the entanglement entropy stems from 
the fact that it quantifies the number of entangled bits 
(or ebits) within the state [11119131] . This means that if 
one has n copies of an entangled pure state, with entan- 
glement entropy S, then one can asymptotically convert 
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this to approximately nS maximally entangled states. It 
has been shown that the maximum amount of shared in- 
formation that two parties can extract from an entangled 
pure state is given by entanglement entropy of their state 
|12I13I14|T5] . Our results thus enable us to determine the 
maximum amount of shared bits per photon pair that two 
parties can extract from SPDC states. Using our results, 
one can also see how changing the parameters of the pump 
or the crystal affects the amount of shared information. 

A straightforward calculation shows that the entangle- 
ment entropy for the SPDC state can be approximated by 
a logarithmic relation; 



7r 



S{K) = lim iJ„ 

a— >1 



l+l0g2(i^). 



(25) 



Such relation is an approximation that holds well for large 
enough values of K . In fact, a Taylor expansion of the non 
approximated Von Neumann entropy for large K yields 



S{K) 




2 + log2if 



1 



/sTlogS 



O 



(26) 



Without considering such power expansion, it is not imme- 
diately obvious how much the relation between Schmidt 
number and Von Neumann entropy differs from a purely 
logarithmic approximation. As an example of the failure 
of the approximation (251, note that in the regime where 
the state is not entangled, i.e. if there is only one joint 
mode (and no K — 1) the amount of quantum correla- 
tion in the state has to be (dashed line in Fig. [T]) and 
not 1 (solid line in Fig. [T]) . Common experimental condi- 
tions where one seeks high dimensional entanglement are 



in the range of &(T <C 1, where the relation (25) is accu- 



rate. However, experiments with a tightly focussed pump, 
or a long crystal, may fall in the region closer to 6cr ^ 1, 
where it fails. A plot of the Von Neuman entropy and of 
its approximation are given in figure [l] 

An interesting feature of figure [T] is that the results are 
symmetric under the substitution ba — >■ 1/ba. This can be 
explained by reference to equation ([2^) which describes 
the correlations in the two conjugate planes (q'j -I- q'^,) 
and (q'j — q'^). Interchanging ba and 1/ba corresponds 
to 'squeezing' in one plane but 'expanding' in the other: 
the product of the two is constant. Physically this means 
that an experiment with a crystal of length L and Rayleigh 
range z^. is equivalent to an experiment with crystal length 
azr and Rayleigh range of L/a. We expect this symmetry 
to be no longer exact outside the approximation of the 
gaussian phase matching function. 

As we already mentioned, the value of the shared bits 
reaches its minimum for ba — 1. In this regime the state is 
separable, i.e. not entangled. In fact, for such value of ba, 
the Schmidt number K is also 1, which means that the 
state can be written using only one Schmidt mode, which 
is separable by definition. In the rest of the paper we will 
assume a large enough Schmidt number to safely use the 
definition (251 of the entanglement entropy. 



The logarithmic dependence of S{K), for large K, has 
two important consequences. Both are due to the small 



5 4 
^ 3 



1 2 3 4 5 6 7 

bo- 
Fig. 1. A plot of the shared bits of information per photon 
pair, given in eq. (25), plotted against ba = y/ L/2zr (solid 



line) and of the non approximated Von Neumann entropy 
(dashed line). Note that they match for large values of K, 
that is for small or very large values of ba. 



value of the slope of the logarithm curve for large values 
of the argument. The first is that if we can prepare an 
SPDC state with a large number of modes, and thus large 
Schmidt number, any further increase in the number of 
modes will provide only a modest increase in both the 
entanglement and the number of shared bits that one can 
extract. The second consequence is concerned with the 
non-ideal detection of the entangled state and is discussed 
in the following subsection. 

It is interesting to compare this result with the work 
|32| . in particular the isotropic case, but also [33 , where a 
relation between spectral entanglement and a control pa- 
rameter similar to ( 19 ) is found. In particular we note that 



in their assessment of the spectral-temporal part of the 
down-converted state, the entanglement strength has its 
minimum at a value considerably larger than zero, mean- 
ing that spectral entanglement is never small. In our case, 
however, the state is spatially separable when ba = 1. In 
both cases, the control parameters depend upon the char- 
acteristics of the crystal and of the pump beam and this 
could have an implication on the extent of hyperentangle- 
ment. 



4.4 Non-ideal detection 

There can be many sources of non-ideal detection. These 
can range from defects in the measuring apparatus (which 
give rise to cross-talk between channels, accidental coin- 
cidences, dark counts, etc.) to turbulence (that can af- 
fect the propagation of the states), to non-ideal choices 
of the optical elements in the setup (which determine a 
mismatch between the Schmidt modes and the detection 
modes and therefore impairs the ability to detect high 
order modes). We concentrate, in particular, on cases in 
which the number of modes that a measurement appa- 
ratus can detect is less than the number of modes that 
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the source is producing. This type of experimental inac- 
curacy is fundamentally different from the others, which 
have been studied for instance in 34 35J, the difference 
being that in the previous work it was always assumed 
that a detector could have access to all the modes that 
are produced by the source. Detection of entangled states 
by projection onto modes that match the Schmidt modes, 
such as the ones given by eqs. ([l2| and ( 15 ), will yield the 
maximum amount of shared bits. However, if the detection 
basis does not exactly match the Schmidt basis, the effec- 
tive number of Schmidt modes that are measured, Kcs^ 
will be less than K. Provided the state is highly entan- 
gled (i.e. large K), the logarithmic relationship between 
the number of shared bits and the Schmidt number, given 
in equation ( 25 ) , means that even if the fraction of en- 
tangled modes that are detected, rj — Kcs/K, is small, 
this need not be too detrimental to the fraction of shared 
bits S(riK) / S{K) that one can extract. This result may 
seem counter-intuitive, however, the key point to note is 
that information is measured by entropy, not by the num- 
ber of modes [25l27j . For example, the number of differ- 
ent messages that one could encode using 4 modes can 
be described using 2 binary digits, while 8 modes require 
3 binary digits, i.e. the information increases by one bit 
every time the number of different messages doubles. 

To illustrate this idea consider the following exam- 
ple. Suppose that one can generate a state with a large 
Schmidt number, K, but that imperfections in the de- 
tection of the modes means that the number of effective 
modes that can be accessed is only K/2, i.e. 77 = 1/2. The 
number of shared bits will then be S{K/2) = S{K) - 1: 
every time K is halved, one shared bit is lost. As we show 
in figure [2] the reduction in the entropy will be negligi- 
ble for large enough values of K, i.e. for small (or large) 
enough values of ba (Fig. [3]). 



If one is interested in determining the experimental 
parameters needed to retain a certain amount of shared 
bits, it is useful to recast Fig. [2] in terms of ba: 



25%los 




0.05 



0.10 



0.15 



0.20 



bo- 



Fig. 3. A plot of the fraction of shared bits per pho- 
ton pair as a function of ba for three types of measure- 
ments that yield different amounts of loss of joint modes. 
The solid lines are derived from the logarithmic approx- 
imation (251, the dashed lines are derived from the non- 



approximated Von Neumann entropy. 



To give some realistic numbers, even if the detection 
basis allows only half of the modes to be detected, it is 
still possible to retain more than 90% of the shared bits, 
if states with more than ~ 500 entangled modes are used, 
which corresponds to ba ^ 0.02. Values of hundreds of 
entangled modes can be produced within the limitations 
of experimental equipment [36j . 



5 Conclusions 




500 



Fig. 2. A plot of the fraction of shared bits per pho- 
ton pair as a function of K for three types of measure- 
ments that yield different amounts of loss of joint modes. 
The solid lines are derived from the logarithmic approx- 
imation ( 25 1 , the dashed lines are derived from the non- 



approximated Von Neumann entropy. 



We have derived analytical Schmidt decompositions for 
the biphoton state produced using SPDC in both carte- 
sian and polar coordinates for cases when the sine phase 
matching term can be approximated as a gaussian. The re- 
sultant modes exhibit either the orthogonal characteristics 
of Hermite-Gauss modes of width F or the angular and 
radial characteristics of Laguerre- Gauss modes of width 
/^/■\/2, respectively, and can be shown to be equivalent 
using the relation between HG and LG modes. An en- 
tropic analysis of these different derivations confirms that 
the strength of the spatial entanglement is independent of 
the choice of coordinates and gives the freedom to choose 
the Schmidt basis that is most appropriate for given ex- 
perimental conditions. 

In order to analyse the entanglement we calculated 
the Schmidt number and the Renyi and von Neumann en- 
tropies for the SPDC biphoton. We showed under what 
conditions it is safe to use a logarithmic relation between 
the number of Schmidt modes and the maximum number 
of bits per photon that one can extract. Such logarithmic 
relation demonstrates that, for highly entangled states, 



Miatto et al.: Cartesian and polar Schmidt bases for SPDC photons 



7 



the information loss due to non-ideal measurements can 
be only a small fraction of the maximum information. We 
finally showed, given any non-ideal measurement, which 
will allow to detect a smaller number of entangled modes 
ifeff, what is the Schimdt number K that is needed in or- 
der to retain a given fraction of shared bits, and from this 
what is the experimental parameter ha that one should 
look for. 



A Cartesian decomposition 

To perform the decomposition in cartesian coordinates we 
require the following mathematical result 



'Y,^^'"h^i^x)K{^y) (27) 



n=0 



for I77I < 1, G > and where 



hniFw) 



(n!2"V^)i/2 



where H„{v) are Hermite polynomials. A proof of this 
formula can be easily obtained with the use of generating 
functions for the Hermite polynomials. 

In order to apply this result to the state (pi) we set 
G ^b^ + and tj = (bV^ - l)/{b^a^ + 1). A bit of 
algebra gives the relations 



M = 







ha -I 


{G + ry2) 




ha + 1 



■,r = 



(28) 



Equation ^ can thus be written in the form 



2a 



2»)gii(jaj_) 

"'■■''h^{rq,)h„,irq,)K{rq,^)K{rqs^) 



(29) 



Let h„in{x,y) — h„i{rx)hn{ry). Using the properties of 
Hermite polynomials one can verify that these functions 
form a complete orthonormal set for £^(IR^) and the result 
m eq. ^ follows. 



B Polar decomposition 

Unlike for the cartesian decomposition which was per- 
formed in one step because cartesian orthogonal degrees 
of freedom play the same role, for the polar decomposi- 
tion it is necessary to separate the angular variables and 
the radial variables in a different way. The angular vari- 
ables will be separated with the Fourier transform, while 
the radial variables will be separated with a variation of 
formula (27 1. 

As thelirst step, we can rewrite the wave function ([2| 
in polar coordinates and show that it is a function of the 



ference of the angular variables. This fact enforces the 
conservation of 0AM and allows to write it as a sum over 
the Fourier components of the difference of the angular 
variables: 



ip — J\f exp 



1 



- {Pi +pI- '^PtPs C0s((/3, - LPs)) 



(30) 
(31) 



where the sum runs over all integers. The Fourier compo- 
nents are easily found: 

v//^F,(p„p,) = AAe-(''^ + ^)(''-+''^)/|,| [2 [h' - ^) p,p, 

(32) 

Where is the order modified Bessel function of 
the first kind. 

The next step is to decompose each angular eigenfunc- 
tion into a radial superposition of orthogonal modes. The 
mathematical result needed to proceed is: 

(33) 



p=0 



a proof of which can be found in 
tions are given by 



Here the rl^-* func- 



rW(x) 



2p! 



,e-^:rl^lifl)(x2) (34) 



where Lp are generalized Laguerre polynomials. 



We apply the formula ( 33 ) to the functions in eq. ( 32 ) 



s_L j A bit of algebra yields the correct value of the parameter 



and the correct scaling, of the r^^"* [Fp) functions: 

'1-6(7 



1+6(7 



;r = 



(35) 



Notice that the values are analogous to the cartesian case. 

Applying these results and normalizing the radial modes, 
the result in eq. ( 14 ) follows. 
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